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Hybrid chiral condensate in the external magnetic field 
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We study the phase diagram of the Nambu-Jona-Lasinio model in the external magnetic field 
within the mean-field approximation, taking into account the inhomogeneous chiral condensate. It 
is shown that there appears a new type of the chiral condensate, endowed with two features of real 
kink crystal and dual chiral density wave, in the magnetic field. We also show that there are first 
order phase transitions between different inhomogeneous phases in the presence of magnetic field. 


I. INTRODUCTION 

In the last decade, the possible appearance of the inhomogeneous chiral phase in the QCD phase diagram has been 
studied jlj-123}. where the quark condensate is spatially modulated. For the analysis of the inhomogeneous chiral 
phase, the chiral order parameter, M(x) = —2 G + i (tW7 5 T 3 ?/>)], has been used. Using the effective models of 

QCD including the Nambu-Jona-Lasinio model or the Schwinger-Dyson approach, there appears the inhomogeneous 
chiral phase in the vicinity of the chiral transition and its critical point is changed to be the Lifshitz point [6|. The 
dual chiral density wave (DCDW) or the real kink crystal (RKC) has been often used as a typical condensate with one 
dimensional spatial modulation. DCDW is a plane wave configuration, M (x) = me lqz , while RKC is a multi-soliton 

configuration, M(x) = (3^7=) ^ > without the phase degree of freedom. In Ref. Q , it has been shown that 

one or two dimensional modulation can be embedded in 1+3 dimensions by using the Lorentz boost. The general 
solutions have been obtained by using the NJL2 model in 1+1 dimensions [Tj, which is called complex kink crystal. 
Higher dimensional modulations has been also considered in some studies pJlElEEIIl- 

In Refs. mi a it has been shown that the inhomogeneous chiral phase can appear at low temperature and 
moderate density region as an intermediate phase prior to the chiral transition. Such inhomogeneous phase has been 
also studied in condensed matter physics, e.g. spin or charge density wave and the FFLO state of the superconductivity 
[24|-[27}. The appearance of the inhomogeneous chiral phase in the QCD diagram has been extensively studied by 
using the various approaches, but there are few works about the external field, isospin asymmetry, and current quark 
mass. These effects should be very important in realistic situations, especially for compact stars. The effect of the 
current quark mass has been studied in some papers [Hi. fiul. [l7l] and change of the phase diagram has been figured out. 

In this paper we consider the inhomogeneous chiral phase in the presence of the external magnetic held to figure 
out some magnetic properties. The effect of the magnetic Held is theoretically and phenomenologically interesting 
and important, since quark matter is put into the strong magnetic Held in compact stars or in heavy-ion collision 
process [28} . QCD in the external magnetic field has recently attracted great attention, and it has been shown that 
the magnetic field gives rise to various phenomena such as chiral magnetic effect [29|, magnetic catalysis [30]-[32|, 
magnetic inhibition [33l - [36} . Lattice QCD simulations have been also performed to study the properties of the QCD 
vacuum in response to the magnetic field jjj ; the effect of the magnetic field on the chiral transition or deconfinement 
has been studied at chemical potential qi = 0. 

The property of the inhomogeneous chiral condensate in the magnetic field has been first studied by Frolov et al. 
[2l[ . They have found that the DCDW phase develops in a wide density region at T = 0 under the magnetic field, 
and that some peculiar behaviours of the amplitude and of the wavevector can be seen due to the de Haas-van Alphen 
effect [13, HI}. However, they did not take into account the possibility of the RKC suggested to be favored in the 
absence of the magnetic field §. 

In this paper, we study the QCD phase diagram in the Nambu-Jona-Lasinio (NJL) model [39l - l4i[ in the magnetic 
field, taking into account both of the condensates. We introduce the a new type of the condensate called hybrid chiral 
condensate (HCC) , M (x) = 1 2 ^L sn ( 1 2 ) ^ 7 , vj e l?z , which smoothly connects both DCDW and RKC by changing 

the modulus v or the wavevector q , and demonstrate that the magnetic field favors the phase modulation: it is found 
through the analysis of the thermodynamic potential that the wavevector q takes a nonzero value in the presence 
of the magnetic field, and thus DCDW and RKC coexist as HCC in the weak magnetic field at moderate densities. 
We shall see that the phase degree of freedom in HCC plays an important role in the presence of the magnetic field. 
The energy spectrum of the quark field becomes asymmetric in the presence of the magnetic field, which gives rise 
to anomalous quark number [2j. Such spectral asymmetry is closely connected with chiral anomaly and moves the 
Lifshitz point to zero chemical potential p, = 0 . 

We consider only the case of isospin symmetric matter (/ i u = fiy) in the chiral limit (m c = 0) for simplicity. In Sec. 
II we briefly summarize the general framework to deal with the inhomogeneous chiral phases in the presence of the 
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magnetic field. We introduce HCC in Sec. III. Spectral asymmetry in the HCC phase and some topological features 
are also discussed there. The phase diagram is presented in Sec.IV in the B — p plane. Sec.V is devoted to concluding 
remarks. The proper-time regularization method is given Appendix A. Some details about spectral asymmetry are 
presented in Appendix B and C. An expansion of the thermodynamic potential with respect to B is given in Appendix 

D . 


II. MODEL AND ENERGY SPECTRUM 

Here we briefly summarize the general framework to get the quark energy spectrum in the inhomogeneous chiral 
phase in the presence of the magnetic field. First, we consider the case of Nf = N c = 1 for simplicity. The case of 
three colors and two flavors is also calculated in the same way. Taking the magnetic field B along , e.g., 2 axis, the 
Lagrangian reads 


L 






1 + 7 " 


M - 


1 -7 5 


M* 


il> 


\M \ 2 

4 G ’ 


( 2 . 1 ) 


within the mean-field approximation, where i/> is 4-dimensional spinor, M is the order parameter of chiral transition, 
M(x) = — 2G [(ipip) + * ( ^7 5r3 '*/’)], ,D^ is covariant derivative D^ = <9^ — ieA y . We consider one dimensional 
modulation along the 2 : axis as well, M(x) = M(z). It is assumed that magnetic field is uniform and parallel to 
modulation of the order parameter. We shall see later that this orientation should be most favorable due to the 
topological aspects [1]. We combined Nickel’s method [H| and Frolov’s method [21] for obtaining the energy spectrum. 
We choose the Landau gauge, = (0, A), A = (0, xB, 0), and assume eB > 0. The Hamiltonian then renders 

H d = a • II + 7°M —+ ^M* —, (2-2) 

where II is kinetic momentum, = —idi + eAi. The Hamiltonian satisfies the commutation relation, [Hjj,—id y \ = 
[Hd, (a_L • n ± ) 2 ] = 0, where aj_ = (a x ,a y , 0) and II j_ = (n^,, n y , 0) . The eigenspinor of — id y and (aj_ ■ HjJ 2 can be 
written as 
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where 77 = xyfeB + k/VeB, u n (rj) is the Hermite function T2] which satisfies (l + iv) u n ( V ) = V2^u n - 1 ( 77 ) and 

— irj'j 1 ( 77 ) = —\[2nu n (y), and n = 0,1, 2,.... denotes the discrete Landau levels. Using this eigenspinor, the 
Hartree-Fock equation = £tjj is reduced to 
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for n = 0. The latter equation (12.51) resembles the Bogoliubov-de Gennes (BdG) equation in 1+1 dimensions, while 
the former equation (12.41) is the same form as the one without magnetic field. Thus the energy spectrum for these 
equations can be obtained once the corresponding one is given in the absence of the magnetic field. The energy 

spectrum in the case of B = 0 is simply written as £ = A±^/1 + p^_/ A± with the perpendicular component of the 
momentum, p± [5;j , where A+ is the eigenenergy of the 1+1 dimensional Hartree-Fock equation Q, 
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and A_ is defined as the eigenenergy for the complex conjugate transformation : M(x) —> M(x)*. Since Eq. (12.411 has 
a similar form to the usual Dirac equation with momentum p± , the eigenvalue can be simply obtained by replacing 
pu by \j2eBn. Thus we obtain the energy spectrum for Eqs. CUD and (g3D ; 


, = / A C V /1 + ^ n = l,2.... 

' 1 A C =+ 


where is asymmetric respect to A<^ = 0 if complex conjugate symmetry is broken, M(x) ^ M(x)*@]. 

These results can be easily generalized to the case of Nf = 2. Assuming that the ground state is the charge 
eigenstate, 


M = —2 G[(ipip) + *{^i 7 5 r 3 ^}] (2.8) 

= {ipi^ 5 T 2 ij)) = 0, (2.9) 


the NJL Lagrangian with three colors and two flavors is written as 
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which is flavor diagonal, so that we can calculate the energy spectrum for each flavor. 
Thermodynamic potential is now written as 
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where p^( A) is the density of states, p^( A) = J]a c ^(A — A^). 


III. HYBRID CHIRAL CONDENSATE IN THE EXTERNAL MAGNETIC FIELD 


We introduce the hybrid chiral condensate (HCC) which has the properties of both of DCDW and RKC, 
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(3.1) 


and is characterized by three parameters; m, q , v. It is reduced to the pure DCDW in one limit, v —> 1, while to the 
pure RKC in the other limit, q —> 0. Thus HCC is the minimum configuration which includes both of DCDW and 
RKC. Note that HCC is simply given by the product of the two types of the condensate, but it satisfies the BdG 
equation within the NJL 2 model. In Refs. [2,0], Basar et al. have found the general form of the condensate in 1 + 1 
dimensions, 


M(z) = — m e iqz A ^ ^ exp [imA z (-i C(i0/2) + ivs(id/2)) + i 9^/2) (3.2) 

a(mAz + 11s. )<j(i 6 /2) 

characterized by four parameters; m,q,u,9, where A = A(9,v) = — 2isc{i9/A\ v)nd(i9/A- : v), and a and ( are Weier- 
strass sigma and zeta functions, and 773 = C(iK / ). When 9 = 2K(i/), this condensate becomes HCC. It can be easily 
seen that the energy spectrum in the HCC phase A^ is uniformly shifted by £g/2 from the RKC one, Af —> — C,q/2. 
Accordingly, the density of states is given as 


Pc(A) = Prkc(A — C<?/2) 
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by using the density of states for RKC, prkc(A) given in 7}, where c = (1 — v — 2E(zz)/K(z/))/(l + \fv) 2 , v' = 
(1 - v^) 2 /(l + \fy) 2 and E(V), K(i/) are the complete elliptic integrals. 



FIG. 1: Profile of HCC on the base manifold, which is the direct product of the horizontal 2 axis and the vertical chiral circle 
given by the scalar and pseudo-scalar condensates. 


Putting the density of states (1331) in Eg. (12. Ill) , we have the thermodynamic potential, which is decomposed into 
the vacuum contribution, the medium contribution and the thermal contribution: 
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Here the vacuum contribution H vac is divergent, so that we use the proper time regularization. 

Hvac = N c lH^/2 ^J/ 5 coth ( T \ e f B \) J dAp c (A)exp (-rA 2 ) 


(3.4) 

(3.5) 

(3.6) 

(3.7) 


IV. SPECTRAL ASYMMETRY AND ANOMALOUS QUARK NUMBER DENSITY 


The fermion number is given by 
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where the first term is the fermion number from spectral asymmetry j43] characterized by the Atiyah-Patodi-Singer 
^-invariant [idj, Qjj which is written as 
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in our case. Here, we have used the fact that spectral asymmetry appears only in the spectrum of the lowest Landau 
level (n = 0), and the higher Landau levels (n yf 0) have no contribution to spectral asymmetry. Note that the 
spectrum becomes symmetric without magnetic field. The second term in Ea. (|4.1l) counts the number of states for the 
given /i with the Fermi-Dirac distribution functions and is usual number density: all the Landau levels contribute to 
this term. Considering quark number density is related to the thermodynamic potential through the thermodynamic 
relation (13.41) : Nb/V = —dQ/dfi, we can explicitly verify that this derivative of thermodynamic potential and Eci. (14.1l) 
are equivalent by using Eq. (1531) . 

The density of states of the lowest Landau level (LLL) is schematically shown in Fig. [2] When the Hamiltonian 
is symmetric for complex conjugation operation M(x) —> M(x)*, the spectrum is symmetric. So in the case of RKC 
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FIG. 2: The behaviour of p+(A). E\ = —m + qj 2, E 2 = —mi/' + q/ 2, E 3 = mv' + q/ 2, E 4 = m + qj 2. DCDW has the gap 
between A = —m + q/2 and A = m + g/2, and has no mid-gap states. The spectrum of DCDW is asymmetric with respect to 
0. The spectrum of RKC is symmetric and has the mid-gap states. In the case of HCC, spectrum is asymmetric and has the 
mid-gap states. 


or homogeneous condensate, number density becomes usual number density. In the case of HCC or DCDW, the 
contribution of //-invariant is nonzero in the presence of the magnetic field. The //-invariant has been evaluated to 
give 
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for DCDW in the case of — m + q/2 < 0 [§]■ It has a topological origin, and is equal to the expression given by chiral 
anomaly [45| (see Appendix iBl). It is straightforward to evaluate the //-invariant for the case of HCC (Appendix [Cj, 
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for the case of —in + q/2 < 0, where the second term is the contribution from the mid-gap states. In particular, for 
to > q/2 > mv ’, it equals to the number of nodes of HCC, N no ^ ea = ni/[(l + v^)K(^)], independent of q. Note that 
this is the same form as in the DCDW phase, except the number of nodes of the condensates. For the general case, 
the / 7 -invariant can be written as //# = VN c J2f ^T’ [~ f + ^pRe (F(g/2; 1 /') + cF(q/2; v') — E(g/2; z/))] , where 
F(x; 1 /') and E (x; v') are the incomplete elliptic integrals of first and second kind. Note that r]n is reduced to the 
DCDW one (14.31) as q —>■ 0, where the energy spectrum is reduced to the one of RKC and symmetric about zero. 

For low chemical potential, to + q/2 > //, usual number density is zero, and the inhomogeneous phase is forbidden 
by the Lorentz symmetry of the vacuum. On the other hand, spectral asymmetry gives nonzero number density, and 
the appearance of the inhomogeneous phase is allowed in the presence of the magnetic field. The contribution of 
spectral asymmetry is taken in the thermodynamic potential as the term, / it)h/2 . Since this term includes the linear 
term in q and the order parameters are determined by the stationary conditions for the thermodynamical potential, 
q = 0 is never the optimal point. In other words the RKC phase itself does not appear in the QCD phase diagram in 
the magnetic field. 

In Q, spectral asymmetry has been also calculated using the derivative expansion and //# = N c Y/f |e/|B • q/2 has 
been obtained for DCDW. Consequently the q ■ B term should appear in the thermodynamic potential by way of the 
thermodynamic relation, and q is favored to be parallel to B. For HCC, the derivative expansion can’t be directly 
applied because the condensate has nodes and the premise that the amplitude is much larger than the wavevector 
breaks down. However, we can manage to evaluate the //-invariant by separating the small nodal region to find 
//h = |e/|B • q/2 + fV no d es for mi/' < q/2, which suggests q is most favored to be parallel to B as well in the 

HCC phase. 


V. PHASE DIAGRAM 

For obtaining the phase diagram, we numerically search the minima of the thermodynamic potential with respect 
to the order parameters; to, q and 1 / for given values of the magnetic field B and chemical potential //. In this paper, 
we show the phase diagram at zero temperature. We use GA 2 = 6.35 and A = 660MeV which reproduce f^ = 93MeV 
and the constituent quark mass ~ 330MeV in the vacuum. 



















6 


A. RKC phase in the magnetic held 

First, we consider the effect of the magnetic field on the RKC phase. Without the magnetic field, the RKC phase 
is energetically more favorable than the DCDW phase in the framework of the NJL model Q. For RKC, the order 
parameter M{z) is real and there appears no spectral asymmetry. 

The phase diagram of the RKC phase is essentially unchanged, while ,as we shall see later, that of the DCDW phase 
is significantly changed in the presence of the magnetic field. Since the chiral condensate is neutral, one may expect 
that there is little effect of the magnetic field on the inhomogeneous chiral phase. Although it holds for the RKC 
phase, some anomalous effect coming from spectral asymmetry plays an important role in the DCDW phase. Figure 
[3] shows the phase diagram of the RKC phase in the presence of the magnetic field. We can see some oscillation of 
the phase boundary with respect to the magnetic field, which comes from the Landau quantization and related to the 
de Haas-van Alphen effect [371|38|. This oscillation is also observed in the case of the homogeneous chiral condensate 
within the NJL model [46| . 

Figure [I] shows the order parameters as functions of chemical potential in different magnetic fields. At y/eB = 
60MeV, the order parameters behave like continuous functions of the chemical potential, and are very similar to those 
in the absence of the magnetic field. At y/eB = 120MeV, the order parameters exhibit some discontinuous jumps, 
since the thermodynamic potential has some local minima as a function of m and v in this region. 
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FIG. 3: Phase diagram for the RKC phase. 


B. HCC phase in the magnetic field 

Here, we consider the phase diagram by introducing the HCC, which includes both features of DCDW and RKC, 
and use the following approximation instead of fully evaluating the thermodynamic potential. When the magnetic 
field is much weak, the approximation 

fl(B) ~ fl(B = 0) + eM2 (1) (5.1) 

is valid, where the first order correction is written in 

eHH (1) = 1/2 (H LLL,q — H LLL-q ) (5.2) 

(See AppendixjD]). This term is an odd function of q , and vanishes at q = 0, so that this term does not appear in the 
thermodynamic potential for the RKC phase. Only LLL contributes to fiW, while the higher Landau levels (n / 0) 
and LLL contribute to the second and higher order terms. We checked the validity of this expansion for some eB by 
comparing the numerical results with Eq. (15.11) and the full expression Eq. (13.41) . Consequently we have found that 
the phase structure is almost unchanged for y/eB < 0.2A ~ 120MeV. 

The magnetic properties of the DCDW phase in the external magnetic field has been studied by Frolov et al. [2l[ . 
They have shown that the DCDW phase is always favorable than that of the homogeneous chiral condensate in the 
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(a )VeB = 60MeV 



(b )VeB = 120MeV 


FIG. 4: Order parameters as functions of chemical potential in different magnetic field, where k is the wavevector of RKC : 
k = 2m/[(l + y/v) K(i/)] 


presence of B. As is already stated in the previous section, it is shown that the mechanism of superiority of the 
DCDW phase is related to spectral asymmetry of LLL states. 

The phase diagram is shown in Fig[5] In this figure, A denotes the weak DCDW phase, B the HCC phase ,C the 
strong DCDW phase, and D the chiral-restored phase. The triple points appear at (/z, y/eB) ~ (320MeV, HOMeV) 
and (/z, VeB) ~ (350MeV, 30MeV). For the limit eB — > 0 , the weak DCDW phase is reduced to the homogeneously 
chiral-broken phase and the HCC phase to the RKC phase. In the eB ^ 0 region, the order parameter is always finite, 
and there is no homogeneously chiral-broken phase nor the RKC phase if the magnetic field has nonzero strength. 
The phase boundary between the chiral-broken and restored phases moves to higher /z as the magnetic field becomes 
stronger. Thus the magnetic field expands the chiral-broken phase mainly due to the phase degree of freedom. 

Figure [ 6 ] shows the energy surface of the thermodynamic potential in the £ — q plane, where m is set to be the 
optimal value for given £, q with £ = Without the magnetic field, there are two local minima corresponding 

to DCDW and RKC in the energy surface of the thermodynamic potential; RKC is energetically more favored than 
DCDW. The minimum with the homogeneous chiral condensate is smoothly changed to that of RKC, so this phase 
transition is of the second order. There is a competition between RKC and DCDW, and there appears no phase in 
which both of phase and amplitude modulations are large. Once turning on the magnetic field, both minima of RKC 
and DCDW move to the larger q direction, which is caused by spectral asymmetry. 

It has been discussed that the mechanism of emergence of the inhomogeneous chiral phase is the Fermi surface 
nesting [IJ. Complete nesting is realized in 1+1 dimensions, when the wavenumber of condensate k c and the Fermi 
wavenumber Jcf have the relation 2 /cf = k c . In 1+3 dimensions, the nesting effect is incomplete, but its reminiscence 
is left in the DCDW phase pj: k c is large, k c ~ 0(2kp)- The situation is a little changed in the HCC phase. 
Using an approximation M(z) — mcos(kz)e lqz , we can decompose the order parameter into two different component: 
M(z) = in.( e l (<i+ k ) z -p = + e *g- 2 )_ The anc j q_ can no t satisfy the nesting relation simultaneously 

if k ^ 0 and q ^ 0. Thus HCC can not satisfy the nesting relation, and there is no HCC phase without the magnetic 
held. In the magnetic held, spectral asymmetry contributes to the emergence of the HCC phase. Note that non-zero 
value of q is favored by some topological effects in this case, different from the nesting effect. 

The order parameters are shown as functions of chemical potential in Fig[7] Under no magnetic helds q is zero 
everywhere, so that there are the homogeneously chiral-broken phase (£ = 0), the RKC phase (£ 7 ^ 0), and the 
chiral-restored phase (£ = 1) whose result is consistent with Nickel’s result. 

With increasing the magnetic held the DCDW phase (£ = 0,g / 0), the HCC phase (£ 0,q ^ 0) and the 

chiral-restored phase (£ = 1) appear. In the low /i region, q is small compared to other order parameters, so we call 
the DCDW phase in the low region ’’weak DCDW phase.” We call the DCDW phase in the high /z region ’’strong 
DCDW phase,” which is similar to the usual DCDW phase with the wavevector q being sufficiently large. In other 
words we may say that the weak DCDW phase is driven by the topological effect due to spectral asymmetry, while the 
strong DCDW phase by the nesting effect. We can see the second order phase transition between the weak DCDW 
phase to the HCC phase, where the order parameters are continuously changed. At the phase transition between the 
HCC phase to the strong DCDW phase, the order parameters change discontinuously. Thus this phase transition is 













FIG. 5: Phase diagram at T=0. A: Weak DCDW phase B: HCC phase C: Strong DCDW phase D: Chiral-restored phase. When 
eB —> 0 , the weak DCDW phase becomes the homogeneously chiral-broken phase smoothly, and the HCC phase becomes the 
RKC phase. The filled circles represent triple points. 


of the first order. 

As we have already seen in the RKC phase, the order parameters should exhibit the de Haas-van Alphen effect 
(H HIl as a function of the magnetic field. In our results, the corresponding effect cannot be seen, since we have 
discarded the contribution of the higher Landau levels in our approximation. If full order contributions of the magnetic 
field is taken into account, the de Haas-van Alphen effect can appear. Anyway, the oscillation of the order parameter 
should be very small at the weak B region , where one may expect the HCC phase. On the other hand, only the 
DCDW phase appears in the high B region, and the phase diagram becomes the same with Frolov’s results [21]. 


VI. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have studied chiral phase transition in the external magnetic field B, taking account of a new 
type of the inhomogeneous condensate called hybrid chiral condensate (HCC). HCC is then a self-consistent solution 
within the NJL model under the mean-field approximation, and exhibits both features of DCDW and RKC. We 
have seen that the quark energy spectrum becomes asymmetric about zero due to the phase degree of freedom of 
DCDW, and there appear the mid-gap states due to the solitonic property of RKC. Generally speaking, spectral 
asymmetry plays very important roles for appearance of inhomogeneous phase in the magnetic field [2[. In some case, 
the contribution of spectral asymmetry is equivalent to manifestation of chiral anomaly. We have explicitly evaluated 
the Atiyah-Patodi-Singer 77 -invariant for the case of HCC. The energy spectrum has a gap and 7777 is given by the sum 
of the one given by the states above and below the gap and the one given by the mid-gap states; the former does not 
depend on the modulus parameter in HCC to give the same form as in the DCDW case, and the latter is related to 
the number of nodes of HCC. 

We have studied the phase diagram of the inhomogeneous chiral phase in the fi — B plane at T = 0 for two cases. 
First we have explored the pure RKC phase and found that the phase diagram is little affected by the magnetic field, 
while some oscillation due to the de Haas-van Alphen effect can be slightly observed. Next, we have discussed the 
full phase diagram, taking into account HCC. Since spectral asymmetry implies that the phase modulation is always 
favored in the magnetic field, independent of /x, the phase diagram consists of three regions besides the chiral-restored 
phase: the strong DCDW phase, the weak DCDW phase, and the HCC phase. The strong DCDW phase resembles 
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(g )VeB = 0.18A, fi = 0.49A 
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FIG. 6: Energy surface of thermodynamic potential at different eB and different fi 


the pure DCDW phase, but the appearance of the weak DCDW phase is attributed to the anomalous quark number 
caused by spectral asymmetry; actually they disappear as the magnetic field is turned off. Note that pure RKC phase 
never appears once the magnetic field is turned on, and is replaced by the HCC phase. 

In this paper, we have considered only the flavor symmetric case, = fid for simplicity, while u and d quarks 
should have different number in a realistic situation due to different electric charge. Actually cold catalyzed matter 
develops inside neutron stars, where charge neutrality and chemical equilibrium should be established. Thus nonzero 
isospin chemical potential (/q = fi u — fid 7^ 0) is very important in the magnetic field. In Refs. ITsL ITol [23l |. they have 
been studied the phase diagram taking into account isospin chemical potential in the absence of the magnetic field. 

We have also considered cold quark matter (T = 0), while it may be interesting to study how thermal effect modifies 
our findings . Actually it has been recently discussed that the external magnetic field suppresses chiral condensates 
at finite temperature [33, MM- This subject is to be discussed elsewhere. 

We have discussed the phase diagram of the inhomogeneous chiral phase in the chiral limit m c = 0, while it is 
known that the current quark mass m c suppresses the inhomogeneous phase QH EE Eil ■ For the case with the 
magnetic field it has been supposed that the effect of current quark mass defeat the effect of spectral asymmetry and 
that chiral condensate is homogeneous at low fi and low B region 0. 

Finally it should be worth mentioning that HCC may have some implications in the context of the FFLO state 
of superconductivity [2(| [27|. Very recently an evidence of the LO state has been reported, and the Andreev bound 
states are emphasized as a hallmark [I?]] . Usually this subject has been separately discussed for the FF state or the LO 
state. Since it has been shown that there is a duality relation between superconductivity and spontaneous breaking of 
chiral symmetry in 1+1 dimensions [H[, the FFLO state with one dimensional modulation may be similarly treated 
to our subject; the Andreev bound states then correspond to the mid-gap states. Note that HCC satisfies the BdG 
equation and can give the pairing function connecting the LO and FF states smoothly. The coexistence of the FF and 
LO states in the quasi-one dimensional system has been discussed as an appearance of time crystal phase in which 
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FIG. 7: Order parameters as functions of chemical potential in different magnetic field. From left to right: (a) eB = 0 (b) 
%/eB = 060MeV (c) y/eB = 120MeV 


time translation symmetry is spontaneously broken [491 ]. 
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Appendix A: Regularization of fi vac 


Since fl v ac is divergent, we apply the proper time regularization for S2 vac . At p = T = 0, thermodynamic potential 
is written in 

ttvac = “^E^E [ d P* f dXp c (X)ln(£ 2 +pl) 
f ^ nX J J 

For ReA > 0 the equation, 


1 

A x 


1 

(x- 1)! 



c - 1 e- rA , 


(Al) 


holds. 

Thermodynamic potential then becomes 


/ n, C 
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^ EE *^ / dXp ^ X) /. 2 ^ exp [~ T ( A2 + 2 \ e f B \ n )] 

f n, C oo Jl/K z 1 

= N - E J 1/A2 ^otHr\e f B\) J ^ dAp c (A)exp (-rA 2 ) 


Appendix B: Some remarks on spectral asymmetry 


Spectral asymmetry is closely related to axial anomaly in the specific case [|}. In the effective theory of mesons, 
anomalous contribution coming from the Wess-Zumino-Witten (WZW) term is given by j45| 

Swzw = ^ 2 j J d 4 a?/iB-V 7 r° (Bl) 

= a 2 + (n 0 ) 2 

in the presence of magnetic field and chemical potential. Our DCDW configuration may correspond to 

a + iT° = Ue i9X , 


in this context. For this configuration, the WZW term reads 

Swzw = ^2 j ■ q (B2) 

Accordingly the anomalous number density is given by 

U = 4^2 B q ' (B3) 

which is the same form as Eq. (Id.31) . This implies that B || q is always favorable. 


Appendix C: Spectral asymmetry for HCC 

Consider the LLL. The rj invariant is then given by 


i)f? = W c E^E hm f dAp+(A)sign(A) |A| s , 

j 2tt s J —qq 


(Cl) 


with the density of states Eq. (liOD . 

For the case, q/2 < to, for simplicity, the integral (IC1I) is divided into two parts: 


/ Ei pE 3 —q/2 f°°\ r°° 1 

+ / + / dAp+(A)sign(A) |A| _S = to / dx— 

-oo J E2~q/2 J E4 J J 1 


X 2 + c 


* y/(x 2 -l)(x 2 -v') 


( mx+ l) *-( mx -I) 


T -^midgap* 

The second integral in Eq. (|C2[) is the contribution of the mid-gap states, 

r E 3 


(C2) 


N, 


midgap 


IE 2 


<2Ap+(A)sign(A). 


For mu' < q/2 , Ad n i ( i gap is equal to the number of nodes of HCC: N noc j es = m/[(l + y/u) K(i/)]. Using the incomplete 
elliptic integrals, N m id gap is written in 2^ Re (F(g/2; u') + cF(q/2; 1 /) — E(g/2; u')) for mu' > q/2. 

In the following, we consider the first term. Expanding it with respect to q , we have 


— 2 TO 


dx-- 


K yj(x 2 — 1)(X 2 — u') 


s(mx) 


-(»+i )£ 1 I 


+ ^s(s + l)(s + 2)(ma;) (s+3) + 0(q 5 ) 


(C3) 
















12 


Since other terms become zero as s —> 0, we, hereafter, evaluate only the first term, 

— sqm~ 


’- s ' dx— 


1 1 7 T yj ( X 2 — l)(x 2 — v') 
Transforming the integration variable x by t = a: -2 , we have 

—— / dtt s/2_1 (l - t) _1/2 (l - ^t) _1 / 2 (l + ct) 
27r Jq 


(C4) 


-sqm 


= -sqm ( 1 / 2 , s/ 2 , (l + s)/ 2 ;l /') + (regular terms in s), 


(C5) 


in terms of the Gauss hypergeometric function F. Using the relation, T(s/2) = 2/sT(s/2 + 1), and taking the limit 
s —^ 0 , we have 


— —F(l/2, 0,1/2; v') = — — 

7T 7T 


Thus the rj invariant can be given as 




midgap J ? 


(C 6 ) 


(C7) 


for the case, q/2 < m. It is easy to evaluate the 77 -invariant for q/2 > m. 

Figure [5] shows the behavior of the number of the occupied states in LLL, -ZVlll, which consists of the normal 
baryon number density and 77 -invariant. For DCDW, the value of JVlll in the plateau is independent of m. For HCC, 
the value of IVlll in the plateau depends on m and v. When the number of nodes is fixed, the value of JVlll in the 
plateau depends on only q. 




(b)HCC (the number of nodes is 
fixed) 


FIG. 8 : The behavior of N\aa i /V as a function of /i for different conditions. 


Appendix D: Expansion of thermodynamic potential with respect to the magnetic field 

Before summation over the Matsubara frequencies, uji = (2 n + 1 )ttT , the thermodynamic potential can be written 


as 


/dAp C (A)I 


\e f B\ 


2tt 

f l,n t C' 




U Z d Vc(A)(l - b„,o)ln [(Wi -i[i) 2 + X 2 + 2\e f B\n] 

f 1 „ /■ J 


f 

— ^even T ^odd- 


ejE 

27r 


~Z r7Vc E ^t^E [ dX ip+W - P-W} ln u? + ( x ~pY 


(Dl) 
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Only LLL contributes to fl 0 dd that is the first order term of q. 

tteven = - \tN c £ EE ^ f dX Pc (X) (l- In [(w, - i/z) 2 + A 2 + 2|e/B|n] 

/ ^ i.n.C' 7 A ' 

= _ 1 TA ^E E /A»cW (1 - r * e -4<"'-M > « , +’lv*l~] (D2) 

2 ^ l,n,C J V y Jo T 

= "^E^E / dA ^c(A) /' 00 ^coth(-r|e / S|)e-[^-^ 2+A2 ]. (D3) 

/ be J Jo 

Since a;coth(a:) is the even function of x, Eleven contains only even order terms of B. 
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